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FOREWORD 


This  report  presents  the  work  accomplished  under  Contract  Nonr-3866(00) 
for  the  period  V~>  December  1962  to  1J>  June  1963. 

The  research  was  carried  out  in  the  Energy  Conversion  RkD  Laboratories  of 
the  Martin  Marietta  Corporation,  Nuclear  Division. 

The  program  was  carried  out  under  the  auspices  of  the  Advanced  Research 
Project  Agency  and  the  Power  Branch  of  the  Office  of  Naval  iiesearch  and  was 
sponsored  by  Dr.  J.  Huth  of  ARPA  and 'Mr.  J.  A.  Satkowski  of  ONR. 
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I.  INTRODUCTION 


The  objectives  of  this  program  are  two:  (1)  to  de"onstrate  experimentally 
the  principle  of  non-equilibrium  ionization  under  actual  closed  cycle  magneto- 
plasmadynamic  (MPD)  generator  operating  conditions;  and  (2)  to  verify  the  basic 
phenomena  underlying  the  behavior  of  moving  ionized  gases,,  in  the  non-equilibrium 
state,  when  interacting  with  a  magnetic  field.  This  will  be  done  by  correlating 
experimental  results  with,  the  theory  of  non-equilibrium  ionization. 

To  bring  about  these  objectives,  a  three-part  program  is  being  carried  out. 

In  the  first  and  major  part  an  electrically  heated,  closed  loop,  experimental 
MPD  electric  power  generator  has  been  developed  and  fabricated  for  the  purpose 
of  evaluating  the  increase  in  output  power  density  and  electrical  conductivity 
due  to  non-equilibrium  ionization.  It  will  also  be  used  to  evaluate  the  dependence 
of  power  generated  in  the  MPD  duct  upon  the  thermodynamic  state  of  the  working 
fluid,  the  magnetic  field  strength,  tne  electrode  temperature,  and  the  electrical 
load  conditions.  The  working  fluid  will  be  a  noble  gas  (helium  seeded  witb  cesium, 
has  been  selected  for  the  initial  experiments),  and  it  will  operate  at  temperatures 
in  the  MPD  generator  duct  within  the  approximate  range  of  1200  to  2000°K. 

In  support  of  the  experimental  part  of  this  program  and  to  guide  in  the 
selection  of  experimental  conditions,  as  well  as  to  correlate  the  experimental 
results  with  theoiy,  two  study  programs  are  being  carried  out:  one  involving  the 
process  of  non-equilibrium  ionization;  and  the  other  involving  the  study  and 
analysis  of  closed  cycle  MPD  power  generation  under  conditions  compatible  with 
the  utilization  of  non-equilibrium  ionization  in  the  MPD  duct. 

All  three  parts  of  the  program  are  thus  directed  toward  demonstrating  the 

success  of  non-equilibrium  ionization  under  actual  closed-cycle  MPD  power  gen- 

0 

erator  operating  conditions.  It  is  also  aimed  toward  gaining  an  understanding 
of  the  basic  important  factors  underlying  the  behavior  of  ionized  seeded  noble 
gases  in  the  non-equilibrium  state,  and  the  application  of  these  principles  to 
MPD  power  generation  at  gas  temperatures  compatible  with  feasible  nuclear  reactor 
capabilities. 
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II.  RESUME  OF  PROGRESS 

During  the  first  six  months  of  the  program'*'  an  experimental  closed  loop 
MPD  electrical  power  generator  was  designed,  and  its  major  components  were  oro- 
cured.  In  addition,  a  theory  on  KPD  electrical  power  generation  w'th  non-equi¬ 
librium  ionization  (of  seeded  noble  gases)  was  developed,  deeded  noble  gas-filled 
discharge  tubes  were  made  and  tests  were  run  to  correlate  experimental  data  with 
the  developed  theory. 

During  the  second  six  months  of  the  program  the  fabrication  of  the  various 
components  of  the  loop :  namely,  the  helium  plenum,  the  helium  flow  meter,  the 
helium  preheater,  the  helium  main  heater  assembly  complete  with  its  six  ceramic- 
to-metal  sealed  main  current  leads,  the  nozzle,  the  MPD  duct  assembly  with  its 
seven  pairs  of  electrod.es,  pressure  and  temperature  probes,  the  main  cooler,  the 
cesium  separator,  the  after  cooler,  the  surge  tank,  and  all  the  interconnecting 
piping  and  valves,  as  well  as  the  instrumentation  and  control  panels,  have  been 
fabricated. 

The  MPD  generator  loop  installation  was  completed,  using  at  first  a  low  tem¬ 
perature  duct  assembly  capable  of  simulated  operation  at  temperatures  in  the 
vicinity  of  800°K.  The  purpose  of  this  low  temperature  duct  assembly  is  to  com¬ 
plete  the  loop  assembly  so  that  initial  leak  checking  and  the  outgassing' of  the 
entire  loop  (as  well  as  the  try  out  of  instrumentation  and  controls)  can  proceed 
without  any  danger  of  contaminating  the  high  temperature  main  heater-nozzle-gen¬ 
erator  duct-diffuser-exhaust  assembly  during  the  process.  The  low  temperature 
assembly  consists  of  a  gas  flow  resistance  section  used  to  simulate  the  resistance 
to  gas  flow  of  the  main  heater,  a  low  temperature  nozzle  section,  a  low  temperature 
duct  section  lined  with  zirconia  insulation  (dimensions:  0.375”  x  1.125"  x  2.625") 
with  four  pairs  of  Mo  electrodes,  a  diffuser  section,  and  an  exhaust  pipe  section. 

*  Semi-Annual  Technical  Summary  Report  (June  15  -  December  15,  1962),  MID -29 39, 

On  Research  Program  on  Closed  Cycle  MPD  Electrical  Power  Generation  with  Non- 
Equilibrium  Ionization. 
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The  loop  assembly,  using  this  low  temperature  duct,  is  now  undergoing 
evacuation  and  leak  checking.  After  location  and  repairing  of  leaks,  certain 
components  of  the  loop  will  be  outgassed  until  a  minimum  evacuation  pressure  in 
the  loop  is  attained.  The  loop  will  then  be  charged  with  the  noble  gas,  the 
impurities  monitored  at  gas  temperatures  up  to  800°K  in  the  duct,  the  loop  purged 
of  the  gas  and  recharged  until  the  level  of  impurities  is  at  a  minimum.  The  loop 
will  then  be  operated  in  try-out  runs  at  a  gas  temperature  in  the  duct  of  800°K 
to  test  the  operation  of  the  loop  controls  and  determine  the  response  of  the  in¬ 
strumentation  to  variations  in  operating  conditions.  The  operation  of  the  cesium 
loop  will  also  be  checked  out  by  circulating  the  cesium  in  bypass  piping  without 
mixing  it  with  the  helium  gas.  After  the  satisfactory  completion  of  leak  checking 
and  the  try-out  phase,  the  low  temperature  duct  assembly  will  be  removed;  and  the 
high  temperature  generator  duct  assembly  will  be  installed,  evacuated,  outgassed, 
purged,  and  charged  with  highest  purity  helium  gas.  The  apparatus  will  then  be 
ready  for  power  generation  during  non-equilibrium  ionization  runs  using  stepwise 
raising  of  gas  temperature  and  cesium  seed  injection. 

The  high  temperature  generator  duct  assembly  consists  of  the  tantalum-10/^ 
tungsten  main  heater  assembly  (with  six  ceramic-to-metal  sealed  main  current  leads), 
the  high  temperature  nozzle,  the  MFD  generator  duct  assembly  (with  alumina. wall 
lining  and  seven  pairs  of  iridium  electrodes),  pressure  and  temperature  probes 
(with  seven  pairs  of  ceramic-to-metal  sealed  output  power  terminals),  diffuser 
section,  and  exhaust  section.  .  These  components  are  all  protected  by  three  series 
of  stainless  steel  cover  boxes  which  will  be  charged  with  a  noble  gas  to  prevent 
oxidation  of  the  refractory  metal  components  while  operating  at  high  temperatures. 

During  the  past  six  month  period,  the  heater  design  was  improved  to  ensure 
against  any  high  temperature  reaction  between  the  tantalum-lC$  tungsten  heater 
elements  and  any  ziroonia  insulation  in  their  proximity.  In  addition,  certain 

design  improvements  have  been  made  and  carried  out  1)  to  ensure  that  the  cover 
boxes  and  current  and  probe  lead-throughs  shall  remain  leak-tight,  and  2 )  to  prevent 
corrosion  or  a  reaction  with  cesium  should  leaks  between  the  gas  flow  and  cover 
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box  gas  regions  occur.  These  design  improvements  and  added  precautions  were 
deemed  necessary  to  allow  longer  trouble-free  runs  with  cesium.  The  major 
assembly  -  consisting  of  the  heater  assembly,  nozzle,  duct  assembly,  diffuser, 
and  exhaust  sections,  along  with  the  protective  cover  box  assemblies  -  is  now 
undergoing  final  leak  testing.  It  will  be  readied  for  installation  after  we 
have  leak-checked  the  rest  of  the  loop  and  have  conducted  satisfactory  trial 
runs  (referred  to  in  previous  paragraphs)  with  the  low  temperature  duct  assembly. 

In  the  second  part  of  the  program  (involving  ionization  studies)  the  theory 
of  "MPD  Electrical  Power  Generation  With  Non-Eowilibrium  Ionization"  was  applied 
in  "A  Comparison  Between  Test  and  Calculations  of  Non-Equilibrium  Conductivities 
and  the  Corresponding  Internal  Electric  Fields  Required"  .  Excellent  agreement 
was  obtained  between  tests  and  calculations  based  on  the  theory  of  non-equilibrium 
ionization  developed  during  the  first  quarter  period  of  the  contract.  In  this 
theory  the  electron  density  is  derived  from  an  ion  balance  equation,  and  the 
electric  field  from  an  energy  balance  equation.  Both  theoretical  calculations 
and  experimental  results  obtained  with  a  diode  tube  (containing  helium  seeded 

p). 

with  about  0.5%  cesium;  helium  density  of  1.61  x  10  atoms  per  cubic  meter)  show 

that  electrical  conductivities  from  10  to  about  15 0  mho/meter  are  obtained  at 

current  densities  in  the  range  of  1  x  lcA  to  20  x  10U  amperes  per  sq  meter  and 

2  2 

with  electric  fields  in  the  range  of  ?  x  10  to  12  x  10  volts  per  meter,  re¬ 
spectively.  These  values  are  of  vital  interest;  and  if  they  are  realizable  in 
MPD  electrical  power  generation  with  non-equilibrium  ionization,  it  would  lead 
to  output  power  densities  on  the  order  of  tens  to  hundreds  of  megawatts  per  cubic 
meter  of  the  generator  duct,  (For  a  more  detailed  report  on  this  investigation 
see  chapter  IV  of  this  report.) 

See  Chapter  III  of  the  first  semi-annual  report,  MND-2939* 

A  ten  minute  paper  presented  before  the  4th  Symposium  on  the  Engineering  Aspects 
of  Magnetohydrodynamics  in  Berkeley,  California  on  April  11,  1963. 
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In  the  third  part  of  the  program,  involving  the  study  and  analysis  of 
closed  cycle  MFD  power  generation,  the  equations  governing  behavior  of  the 
ionized  gas  in  the  MFC  duct  have  been  solved  for  conditions  compatible  with 
the  utilization  of  non-equilibrium  ionization  in  the  MFD  duct.  These  solutions 
are  presented  in  chapter  V  of  this  report  for  the  case  of  a  Faraday  generator 
with  segmented  electrodes,  operating  with  constant  velocity  or  constant  area 
duct  and  with  constant  electrical  conductivity  along  the  generator  duct  length. 

For  compatibility  with  the  principles  of  non-equilibrium  ionization , the  further 
condition  is  imposed:  that,  in  order  to  maintain  a  certain  conductivity  under 
the  non-equilibrium  state,  it  is  necessary  to  maintain  constant  the  current  density 
between  each  electrode  pair  and  the  corresponding  Internal  electric  field  along 
the  generator  duct.  To  our  knowledge,  this  is  the  first  time  the  MPD  generator 
equations  have  been  solved  for  the  constant  area  duct  case  under  conditions  com¬ 
patible  with  the  generation  of  non-equilibrium  ionization  conductivities.  The 
results  are  presented  in  chapter  V  in  normalized  universal  solutions  to  assist 
in  their  application  to  any  examples  of  duct  design,  in  performance  prediction, 
or  in  correlation  with  theory. 
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III.  EXPERIMENTAL  CLOSED  CYCIE  MPD  ELECTRICAL  POWER  GENERATOR  • 

1.  Introduotion 

The  design  of  the  experimental  closed-cycle  KPD  electrical  power  generator 
was  described  in  detail  in  the  first  semi-annual  progress  report  (KND-2939). 

Copies  of  the  design  drawings  of  the  overall  equipment  layout  of  the  MPD  ex¬ 
periment  (JU35A-6i+0001U),  the  MPD  instrumentation  and  control  block  diagram 
(E43!?A1367000) ,  the  helium  plenum  (EU35-6U00009) ,  the  main  heater-nozzle-duct- 
dif fuser-exhaust  assembly  (Fig.  12  of  the  first  semi-annual  report,  MND-2939), 
the  main  cooler  ( JU35>A6I;OOOOl|) ,  the  cesium  separator  (EU35-6U00006) ,  the  after 
cooler  (EU3561t.00005 )  and  the  cesium  accumulator  (EU35-6)(00015>)  were  included  in 
the  same  report  (KND-2939) .  Pictures  of  the  heater  elements  and  plates  (Figs.  U 
and  11  of  MND-2939),  some  of  the  components  making  up  the  heater  and  duct  assemblies 
(Fig.  13  of  MND-2939) ,  the  compressor  piston  and  piston  rod  (Fig.  lh  of  MND-2939 ) , 
the  electromagnet  (Fig.  16  of  KND-2939)  and  the  building  housing  the  loop  (Fig.  17 
of  MND-2939)  were  also  included  in  the  first  semi-annual  report. 

2.  Progress  This  Period 

Since  the  beginning  of  this  period,  the  helium  plenum,  the  helium  flow  con¬ 
trols  and  helium  flow  meter,  the  helium  preheater,  the  cesium  injector,  the  main 
cooler,  the  cesium  separator,  the  after  cooler,  the  helium  surge  tank  (ahead  of 
the  inlet  to  the  compressor),  the  cesium  storage  tank,  the  cesium  supply  tank, 
the  cesium  pump,  the  cesium  accumulator,  the  heated  oil  bath  tank  (housing  the 
cesium  tanks,  cesium  pump  head  and  accumulator),  the  cesium  flow  controls,  the 
cesium  flow  meter,  and  all  the  interconnecting  piping  and  valves  have  been  fab¬ 
ricated  (or  procured)  and  assembled  in  the  loop.  Compressor  installation  in  the 
loop  was  completed  by  connecting  the  surge  tank  to  the  compressor  inlet,  by  com¬ 
pleting  the  pipe  connection  between  the  stages  of  the  compressor,  and  then  by 
connecting  the  compressor  outlet  to  the  helium  plenum.  . he  variable  speed  com¬ 
pressor  drive  was  aligned  with  and  coupled  to  the  compressor  shaft,  and  its 
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rotation  and  balance  were  tested.  In  addition,  the  instrumentation  and  control 
panel  were  fabricated,  assembled,  and  initiated  in  the  control  room,  located 
next  to  the  loop. 

The  major  subassembly,  consisting  of  the  main  heater-nozzle  -  generator  duct_dif- 

\ 

fuser-exhaust  duct,  together  with  the  six  ceramic-to-metal  sealed  main  heater  current 
leads,  the  fourteen  ceramic-to-metal  sealed  generator  duct  output  power  leads, 
pressure  and  temperature  probes  and  the  three  enclosing  cover  boxes  for  protection 
of  the  refractory  metal  parts,  was  completed  on  May  17,  1963.  This  represented 
a  delay  of  about  four  weeks  from  the  original  schedule;  and  it  was  due  primarily 
to  difficulties  encountered  in  fabrication  of  the  six  ceramic-to-metal  sealed  main 
heater  current  leads  (seventeen  days  delay),  and  to  incorporation  of  design  im¬ 
provements  to  ensure  effective  isolation  of  all  the  tantalum-lOfo  tungsten  parts 
(e.g.  the  heater  elements,  heater  assembly  box,  nozzle  support  structure,  generator 
duct  assembly  box,  etc.)  from  any  zirconia  insulation  blocks  (eleven  days  delay). 

This  design  modification  involved  replacing  the  molybdenum  isolation  sheets  (between 
the  tantalum-10/o  tungsten  material  and  the  zirconia  insulation)  with  tungsten 
isolation  sheets.  This  increases  the  effectiveness  of  the  isolation  layer  as  a 
barrier  to  possible  reaction  between  the  tantalum-10^  tungsten  and  the  zirconia 
materials,  particularly  at  higher  temperatures  (e.g.  above  1300°K).  The  completed 
assembly  was  leak-checked  during  the  week  between  May  17  and  Kay  2U  and  leaks  were 
located  in  some  of  the  welds  in  the  stainless  steel  cover  boxes  and  in  the  tantalum- 
10^  tungsten  exhaust  duct,  as  well  as  in  two  internal  O-ring  seals  between  the  main 
gas  stream  region  and  the  inert  gas  region  inside  the  cover  boxes.  It  was  determined 
that  repairing  the  leaks  in  the  welds  could  be  done  readily  and  without  significant 
delay  while  repairing  the  leaks  in  the  two  internal.  0- rings  would  involve  a  major 
"tear-up"  operation,  accompanied  by  considerable  risk  and  an  extensive  delay  in 
the  schedule.  On  May  27,  therefore,  it  was  decided  to  proceed  with  repair  of  all 
leaks  in  the  stainless  steel  and  tantalum-10$  tungsten  welds.  Moreover,  since  any 
leak  between  the  main  gas  stream  region  and  the  inert  gas  regions  under  the  cover 
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boxes  would  not  be  serious  unless  external  room  air  leaked  into  the  boxes  and 
reacted  with  the  refractory  metals,  it  was  concluded  that  the:  cover  boxes  should 
be  made  leak  tight..  To  ensure  this,  it  was  decided  to  replace  all  brazed  pressure 
and  temperature  probe  or  current  lead-throughs  in  the  cover  box  walls  which  might 
react  with  cesium  (in  case  it  leaked  into  the  cover  boxes  inert  gas  regions  from 
the  main  gas  stream  region),  with  leak-tight  stainless  steel  welded  pressure  and 
temperature  probe  and  current  lead-throughs  using  insulation  materials  which  would 
not  react  with  or  be  corroded  by  cesium.  These  design  improvements  have  been 
carried  out,  and  at  this  writing,  the  major  assembly  (or  the  main  heater-nozsle- 
generator  duct-diffuser-exhaust  duct  together  with  the  six  ceramic-to-metal  sealed 
main  heater  leads,  the  fourteen  ceramic-to-metal  sealed  generator  duct  output  power 
leads,  pressure  and  temperature  probes  and  the  three  enclosing  cover  boxes)  is 
being  rechecked  for  leaks. 

Furthermore,  it  was  decided  to  close  the  loop  temporarily,  using  a  low  tem¬ 
perature  duct  assembly,  in  order  to  proceed  with  the  evacuation,  leak-checking, 
and  outgassing  of  the  remainder  of  the  loop  assembly.  In  addition,  the  use  of 
the  low  temperature  duct  assembly  would  permit  trying  out  the  helium  gas  purging 
process  and  impurity  monitoring  procedure  (using  a  chromatograph)  -  as  well  as 
the  checkout  of  all  the  instrumentation  and  controls  in  trial  runs  at  gas  tern-, 
peratures  (in  the  duct)  of  about  800°K.  This  low  temperature  duct  assembly  is 
made  of  stainless  steel  and  consists  of  a  gas  flow  resistance  section  (to  simulate 
the  resistance  to  gas  flow  through  the  main  heater),  a  low  temperature  nozzle 
section,  a  low  temperature  duct  section  lined  with  zirconia  insulation  (dimensions: 
0.375"  x  1.125"  x  2.625")  and  having  four  pairs  of  molybdenum  electrodes  (0.375"  x 
0.375"  each),  a  diffuser  section,  and  an  exhaust  pipe  section.  The  loop  was  closed 
with  this  low  temperature  duct  assembly  on  June  2h>  and  it  is  now  being  evacuated 
and  leak-checked. 

The  following  figures  (Figs.  1-13)  show  .the  loop  installation  using 
the  low  temperature  duct  assembly,  the  compressor  and  compressor  drive  installation 
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the  instrumentation  and  control  panel,  and  the  major  assembly  of  the  main  heater- 
nozzle-generator  duct-diffuser-exhaust  duct  and  protective  cover  boxes  in  different 
stages  of  completion. 
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FIGURE  CAPTIONS 


Figure  1. 
Figure  2. 

Figure  3. 

Figure  U* 
Figure  5. 

Figure  6. 

Figure  7. 

Figure  8. 

Figure  9. 

Figure  10. 
Figure  11. 
Figure  12. 
Figure  13, 


MPD  Loop  Assembly  -  View  Showing  Low  Temperature  Duct  between  the 
Magnet  Poles 

MPD  Loop  Assembly  -  View  lowing  Exhaust  Pipe,  Cooler,  Separator, 

After  Cooler  and  the  Oil  Tank  Housing  Cesium  Storage  Tanks,  Accumulator 
and  Pump  and  Cooler 

MPD  Loop  Assembly  -  View  Showing  the  Compressor  With  Its  Drive  and 
the  Surge  Tank 

MPD  Loop  Control  Panel 

MPD  Generator  Major  Assembly  Embodying  the  Main  Heater  Leads,  the 
Main  Heater,  the  Nozzle,  the  Generator  Duct,  the  Diffuser,  the  Exhaust 
Duct  and  the  Output  Power  Leads,  and  the  Pressure  and  Temperature 
Probes 

MPD  Generator  Major  Assembly  -  View  Without  Cover  Box  on  Main  Heater 
Lead  Seal  Assemblies 

MPD  Generator  Major  Assembly  -  View  Showing  Output  Power  Leads, 

Pressure  and  Temperature  Probes  Coming  Through  the  Cover  Box 

MPD  Generator  Major  Assembly  -  View  Showing  the  Output  Power  Leads 
and. Refractory  Ceramic-to-Metal  Seals  Inside  the  Cover  Box 

MPD  Generator  Major  .Assembly  -  View  Showing  Method  of  Assembling  the 
Front  Cover  Box  on  the  Main  Heater  Lead  Seal  Assemblies 

Input  Power  Lead  Seal  Assembly 

Heater  Box  and  Element  Assembly 

Generator  Electrodes  and  Insulation  Components 

Take-Off  Seal  and  Flexible  Lead  Assembly 
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IV.  NON-EQUIUBRIUM  IOHIZATIOU  STUDIES 

This  chapter  is  covered  by  the  text  of  the  following  technical  paper: 

"Comparison  Between  Test  and  Calculations  of  Non-Equilibrium  Ionization  Con-  ■ 
ductivities  And  the  Corresponding  Internal  Electric  Fields  Required,"  presented 
at  the  Fourth  Symposium  on  the  Engineering  Aspects  of  Kagnetohydrodynamics  held 
at  Berkeley,  California,  April  11,  1963.  Its  text  represents  pages  11-17  of 
this  report,  and  Figs.  1  to  U  of  the  paper  represent  Figs.  1U  to  17  of  this 
report. 
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COMPARISON  BETWEEN  TEST  AND  CALCULATIONS  OF'  NON-EQUILIBRIUM 
IONIZATION  CONDUCTIVITIES  AND  THE  CORRESPONDING 
INTERNAL  ELECTRIC  FIELDS  REQUIRED 

M.  E.  Talaat 

Martin  Marietta  Corporation,  Baltimore  3,  Maryland 

INTRODUCTION 

An  ionized  gas  is  said  to  be  in  the  non-equilibrium  state  when  the 
approximately  Maxwellian  electrons  are  maintained  at  a  higher  temperature 
than  the  Maxwellian  ions  and  neutrals.  When  such  a  state  exists  in,  for 
example,  a  seeded  noble  gas  the  ions  could  be  primarily  generated  by  ion¬ 
izing  collisions  of  the  hot  plasma  electrons  with  the  seed  atoms  and  the 
gas  temperature  could  be  maintained  at  a  sufficiently  low  value  which  is 
compatible  with  a  nuclear  reactor  heat  source.  The  term  non-equilibrium 
ionization  is  used  to  refer  to  the  ionization  processes  under  these  con¬ 
ditions.  The  realization  of  a  practical  closed  loop  magnetoplasmadynamic 
generator  with  a  reactor  heat  source  depends  on  the  feasibility  of  attaining 
good  electrical  conductivities  by  application  of  the  principles  of  non¬ 
equilibrium  ionization  to  such  a  system. 

Thus,  it  is  important  to  be  able  to  predict  with  close  accuracy  the 
electrical  conductivity  of  an  ionized  gas  in  the  non-equilibrium  state  and 
the  corresponding  electric  field  required  for  maintaining  the  hot  plasma 
electrons  and  in  turn  the  rate  of  generation  of  ions  or  ionization  level 
which  is  necessary  to  yield  that  electrical  conductivity.  Accuracy  in  pre¬ 
dicting  the  electrical  conductivity  is  necessary  because  the  output  power 
density  of  a  magnetoplasmadynamic  generator  is  directly  proportional  to  the 
conductivity,  and  accuracy  in  predicting  the  corresponding  electric  field 
required  is  necessary  because  it  gives  an  indication  of  the  feasibility  of 
attaining  the  desired  conductivity  with  non-equilibrium  ionization  under 

any  given  set  of  operating  conditions. 


In  this  paper  we  present  a  comparison  between  the  calculated  and  measured 
values  of  the  electrical  conductivity  and  the  corresponding  electric  fields 
required  in  a  gas  filled  diode  containing  helium  seeded  with  cesium.  The 
calculations  are  based  on  the  theory  of  non-equilibrium  ionization  which  was 
developed  in  Reference  1,  and  in  which  the  electron  or  ion  density  in  the 
plasma  is  derived  from  an  ion  balance  equation  which  equates  the  rate  of 
generation  of  ions  by  the  hot  electrons  in  the  plasma  to  the  rate  of  loss 
of  ions  by  recombination.  The  electric  field  required  to  maintain  the  hot 
electrons  in  the  plasma  is  derived  from  an  energy  balance  equation  which 
equates  the  rate  of  energy  lost  by  the  electrons  in  elastic  and  inelastic 
collisions  with  the  gas  species  to  the  rate  of  energy  fed  to  the  electrons 
through  the  electric  field. 

ENERGY  BALANCE  EQUATION 

In  the  case  of  non-equilibrium  ionization  discussed  in  this  paper,  the 
rate  of  energy  lost  by  the  plasma  electrons  in  ionizing,  exciting  and  elastic 
collisions  with  the  cesium  seeded  helium  gas  species  is  balanced  by  the  rate 
of  energy  fed  by  the  flow  of  a  current  density,  Ig,  under  the  electric  field, 
£e .  Thus, 

-  Pel  *  Pinel 

If  the  variation  in  the  electrical  conductivity  of  the  plasma  is  negligible, 
the  diffusion  equation  for  electrons  yields  the  following  relation  for  the 
current  density 


and  the  energy  balance  equation  takes  on  the  following  form 
<r£2e  "  Pel  *  Pinel 
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where  for  a  Maxwellian  distribution  of  electrons  in  a  partially  ionized 
cesium  seeded  helium  gas  we  have  (see  Refs.  1  aiid  2) 
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ELECTRIC  FIEID 


Next  we  substitute  the  egressions  for  ^  and  CT  into  the  energy 

balance  Eq.  (3)  and'  we  obtain  the  following  expression  for  the  electric  field, 

£  i  inside  the  plasma  which  is  required  to  maintain  the  flow  of  current  and 

hence  feed  into  the  plasma  the  rate  of  energy  necessary  to  maintain  the 

approximately  Maxwellian  plasma  electrons  at  any  specific  electron  temperature 

E  ,  namely, 

0 

£  >  2.12n  q  E  -/F 
°e  g^g  e  T 

where  f  is  the  equivalent  average  fraction  of  the  electron  energy  which  is 
lost  in  elastic  and  inelastic  collisions  and  is  given  by  (see  Ref.  1) 


ELECTRON  DENSITY 

The  electron  density,  ng,  appearing  in  the  expression  for  the  electric 
field,  E  ,  and  in  the  expression  for  the  electrical  conductivity,  o'  ,  is  a 
function  of  the  electron  temperature  and  is  derived  from  the  generalized  dif¬ 
ferential  expression  of  the  ion  balance  equation  which  equates  the  rate  of 
generation  of  ions,  V  to  the  rate  of  loss  of  ions  by  diffusion  —  V • 

and  by  volume  recombination  o<  n  n  ,  namely, 

rep 
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Now  if  the  conditions  are  such  that  the  diffusion  to  the  walls  is  negli¬ 
gible  by  comparison  to  volume  recombination  and  if  the  gas  temperature,  E  , 

g 

is  low  by  comparison  to  the  electron  temperature,  Eg ,  such  that  the  only  im¬ 
portant  mechanism  of  ionization  (see  Ref.  1)  is  by  collision  of  the  hot 

electrons  with  the  seed  atoms  (viz.,  the  rate  of  generation  of  ions V\ ■  S.ti  n  ) 

,  1  1  s  e 

then  the  ion  balance  equation  yields  the  following  approximate  equation  for  the 
electron  density 


TTT5T  n 


so 


where  the  recombination  coefficient,  °<T,  has  in  most  cases  of  interest  as  the 
dominant  component  the  Thomson  3  body  recombination- coefficient,  c*  ^  ,  namely 
(see  Ref.  1) 

Vos  +  Vp)  # 

Vg  J 


<*  -(2.68)10“21E  '2 


•5V*  (J 


1  +  0.752 


RESULTS 

The  electrical  conductivity  and  the  electric  field  versus  current  density 
curves  for  any  given  case  of  seeded  noble  gas  can  now  be  computed  by  solving 
these  equations  for  the  electron  density,  the  corresponding  electrical  conduc¬ 
tivity,  the  corresponding  electric  field  and  the  corresponding  current  density 

at  several  chosen  values  of  the  electron  temperature,  E  . 

© 

The  results  of  the  calculations  are  presented  in  Figs.  1  and  3  for  the  case 

of  a  helium  gas  density  of  1.61  x  10  4atoms/nr  seeded  with  a  cesium  seed  density 
22  3 

of  0.8l8  x  10  atoms/nr  (corresponding  to  a  liquid  cesium  temperature  of  523°K 
or  0.015  volt).  In  one  set  of  calculations  the  gas  temperature  is  taken  to  be 
the  same  as  the  temperature  of  the  electrodes  in  the  diode  tube,  viz.,  Eg“0.114i 
volt  (corresponding  to  a  temperature  of  the  electrodes  *  l670°K)(see  Fig.  1)  while 
in  another  set  of  calculations  the  gas  temperature  is  taken  to  be  the  same  as  the 
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liquid  cesium  temperature,  viz.,  E  =  .045  volt  (see  Fig.  3).  The  solid  lines 

o 

are  based  on  a  value  for  the  initial  slope  of  the  ionization  cross  section  curve 
for  cesium  derived  from  Von  Engel  k =  0.79  x  10  sq  m/volt  and  with  an  assumed 
equal  value  for  the  initial  slope  of  the  excitation  cross  section  curve  for  cesium. 
The  dashed  curves  are  based  on  a  value  of  the  initial  slope  of  the  ionization  cross 
section  curve  for  cesium  »  10  times  that  used  in  deriving  the  solid  curves,  but 
with  the  initial  slope  of  the  excitation  cross  section  curve  equal  to  that  used 
in  deriving  the  solid  curves. 

The  experimental  points  indicated  in  Fig.  1  were  obtained  from  tests  on  a 
0.0135  meter  long  diode  tube  with  hot  electrodes  and  filled  with  the  cesium-seeded 
helium.  These  points  were  evaluated  from  the  voltage-current  density  curve  shown 
in  Fig.  2  (taken  from  Fig.  27  of  Reference  4)  with  the  electrodes  maintained  at  a- 
temperature  of  1670°K. 

The  experimental  points  indicated  in  Fig.  3  were  obtained  from  tests  on  a 
0.23  meter  long  diode  tube  with  hot  electrodes  and  filled  with  the  cesium-seeded 
helium.  These  points  were  evaluated  from  voltage  and  current  readings  at  a  time 
of  about  25-50  micro  sec  after  the  initiation  of  the  discharge  which  was  long  with 
respect  to  the  time  required  for  collisional  ionization  equilibrium,  but  short  with 
respect'  to  the  time  for  heating  the  gas.  Thus,  all  experimental  points  for  this 
figure  can  be  considered  to  have  been  obtained  at  approximately  the  same  steady- 
state  gas  temperature  of  523°K. 

It  is  seen  that  the  agreement  between  the  theoretical  and  experimental  results 
is  very  good.  The  calculated  values  of  electric  fields  required  to  maintain  a 
certain  current  density  are  generally  higher  than  those  measured  in  the  experiment 
for  current  densities  in  the  range  of  about  1  to  about  20  amps  per  sq  cm.  The 
positive  resistance  characteristic  predicted  by  the  theory  is  also  verified  by 
the  experimental  results.  The  calculated  solid  curve  values  of  the  corresponding 
electrical  conductivities  are  lower  than  those  measured. 
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Since  these  calculations  are  based  on  values  of  electron  densities  which 
are  computed  using  the  ion  balance  equation  rather  than  those  which  would  be 
obtained  from  using  the  electron  temperature  in  Saha's  equation  for  thermal 
ionization,  it  is  of  interest  to  compare  plots  of  the  calculated  electrical 
conductivities  versus  electron  temperature  using  the  two  approaches  in  obtaining 
the  electron  densities.  In  Fig.  I4  the  dotted  curve  was  derived  using  the  electron 
density  from  Saha's  equilibrium  equation  with  the  electron  temperature  inserted. 

The  solid  and  dashed  curves  were  derived  for  E  *  0.0ii5  volt  using  the  electron 
density  as  obtained  from  the  ion  balance  equation  with  the  value  of  the  initial 
slope  of  the  ionization  cross  section  in  the  dashed  curve  taken  to  be  equal  to 
ten  times  the  value  given  in  Von  Engel's  book  on  Ionized  Gases  and  which  was  used 
in  deriving  the  solid  curve.  It  is  seen  that  the  curve  using  the  electron  density 
from  Saha's  equation  with  E  inserted  gives  higher  values  of  electrical  conduc- 
tivities  (sometimes  by  more  than  two  orders  of  magnitude)  than  both  curves  using 
the  electron  density  from  the  ion  balance  equation  with  an  ionization  cross  section 
of  cesium  as  given  by  Von  Engel  or  with  ten  times  that  value. 

CONCLUSION 

In  conclusion,  the  comparison  between  test  and  calculations  based  on  a  theory 
of  non-equilibrium  ionization  (Ref.  1)  in  which  the  electron  density  is  derived 
from  an  ion  balance  equation  and  the  electric  field  is  derived  from  an  energy 
balance  equation  indicates  very  good  agreement  between  experiment  and  theory.  Both 
theoretical  calculations  and  experimental  results  obtained  in  the  diode  tubes  filled 

pi,  o 

with  helium  seeded  with  about  0,5%  cesium  and  helium  density  of  1.61  x  lOatoma/r 
show  that  electrical  conductivities  on  the  order  of  100  mho/m  are  obtained  at  cur¬ 
rent  densities  on  the  order  of  10  x  10^  amp/sq  m  and  with  an  electric  field  on  the 
order  of  1200  volt/m.  These  values  are  of  interest  in  KPD  electrical  power  gen¬ 
eration  with  non-equilibrium  ionization. 
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V.  STUDY  AMD  ANALYSIS  OF  CLOSED-CYCLE  MPD  POWER  GENERATION  UNDER  CONDITIONS 
COMPATIBLE  WITH  NON-EQUILIBRIUM  IONIZATION'  ‘ 

1«  Introduction 

In  the  following  articles  and  subsequent  articles  in  future  reports  under 
this  title  we  will  analyze  the  operation  of  the  closed-cycle  MPD  power  generator 
under  conditions  which  are  compatible  with  attaining  the  electrical  conductivity 
in  the  MPD  duct  by  means  of  non-equilibrium  ionization  techniques  (either  by 
magnetically  induced  electron  heating  or  by  other  auxiliary  ionization  techniques 
such  as  photoionization,  electron  injection  gun,  auxiliary  d-c,  a-c  or  h.f .  dis¬ 
charge,  etc.).  We  will,  first,  in  the  following  sections,  obtain  solutions  of 
the  MPD  power  generator  equations  for  the  following  cases : 

a.  Constant  velocity  segmented  electrodes  Faraday  generator  with  constant 
electrical  conductivity,  (j^,  and  constant  internal  voltage  drop  between 
each  electrode  pair  along  the  duct  (i.e.  constant  current  density  flow 
between  each  electrode  pair  along  the  duct). 

b.  Constant  area  duct  segmented  electrodes;  Faraday  generator  with  constant 
electrical  conductivity,  (J~0,  and  constant  internal  voltage  drop  between 
each  electrode  pair  along  the  duct  (i.e.  constant  current  density  flow 

•between  each  electrode  pair  along  the  duct). 

c.  Constant  area  duct  segmented  electrodes  Faraday  generator  with  constant 
electrical  conductivity,  (T^,  and  constant  electrical  load  resistance 
for  each  electrode  pair  along  the  duct. 

d.  Constant  area  duct  segmented  electrodes  Faraday  generator  with  constant 

electrical  conductivity,  and  constant  output  voltage  along  the  duct. 

The  solutions  (e.g.  power  density,  efficiency  and  duct  length)  are  presented 
here  in  figures  using  universal  dimensionless  factors  which  can  be  used  in  the 
analysis  or  design  of  any  generator  of  the  types  specified  in  cases  (a)  to  (d). 
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In  subsequent  articles  in  future  reports  under  this  title  we  will  examine 
the  solution  of  the  MPD  power  generation  equations  for  the  case  of  segmented 
electrodes  Hall  generator  and  under  conditions  similar  to  those  listed  in  (a) 
to  (d). 

These  solution's  will  be  used  to  evaluate  the  MFD  power  generation  cycle 
in  which  the  reject  heat  from  the  MPD  generator  duct  is  used  to  power  the 
compressor  necessary  for  closing  the  cycle.  In  this  context  we  will  evaluate 
and  assess  the  use  of  the  jet  compressor  as  well  as  the  conventional  mechanical 
compressor.  We  will  also  present  the  basic  equations  governing  the  behavior 
of  the  simple  constant  area  mixing  jet  compressor. 

Examples  of  calculations  will  be  selected  and  numerical  evaluations  of 
the  performance  of  typical  generator  and  its  suitable  compressor  will  be 
carried  out. 
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2.  General  Equations  of  the  MPD  Power  Generator 

The  equations  governing  the  behavior  of  the  KPD  electric  power  generator"*- 
are  Maxwell's  electromagnetic  field  equations,  the  generalized  Ohm's  law  and 
the  gas  dynamic  equations  modified  to  take  account  of  the  interaction  between 
the  motion  and  the  magnetic  field  and  the  energy  equation. 

Maxwell's  Equations; 

Maxwell's  equations  consist  essentially  of  the  following: 

a.  The  generalized  differential  expression  of  Faraday's  law  for  the  curl 
of  the  electric  field,  viz*, 

v  x  6  «-|f-  (l) 

b.  The  generalized  differential  expression  of  Gauss'  law  for  the  divergence 
of  the  magnetic  field,  viz., 

B  — 0  (2) 

c.  The  generalized  differential  expression  of  Biot  and  Savart's  law  for 
the  curl  of  the  magnetic  field,  viz., 

7  x  B  -  /<  •  J  +  /*  Q  (3) 

d.  The  generalized  differential  expression  of  Gauss'  law  for  the  divergence 
of  the  electric  field,  viz., 

(^/€q)  (Poisson's  Equation)  (k) 


1.  Talaat,  M.  E.,  "Magnetohydrodynamic  Electric  Power  Generators,"  Advanced 
Energy  Conversion,  Vol.  1,  pp.  19-35,  Pergamon  Press,  1961. 
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Ohm's  Law 


The  generalized  differential  expression  of  Ohm's  law^  for  a  moving  medium 

includes  beside  the  term  due  to  the  electric  field,  ,  a  term  due  to  the  in- 

J 

duced  electric  field,  (u  -  ~ — )  x  B  ,  resulting  from  the  motion  of  the  medium 

e 

at  velocity,  u,  under  the  influence  of  the  magnetic  field  intensity,  B,  a  term 
due  to  the  transport  of  the  net  charge  density,  A ,  by  the  medium  moving  with 
a  velocity,  u,  and  a  term  due  to  the  electron  pressure  gradient  VPe  “ 


€  *  (u  — — )  x  B  +  -i-  »  p  1  +  P  u 
—  v—  en  '  -  en  v  *e  '  q  — 

e  e  J 


J  -  (T 

Note  that  if  u  =  0  and  B  =  0  then 


(5) 


J  =  C  6  +  v  P  =  en  £  +  A  7  (n  kl  )  »  en  A  6  +  V  (en  £  A  )  (5a) 

—  —  en  e  e  —  e  ee  e  e—  eee 


which  is  the  diffusion  equation  of  electrons  under  no  magnetic  field^. 
When  u  =  0  and  B  f  0  then 

J  -  V  (nekTe)  -  ;i-  J  *  B 

e 

which  is  the  diffusion  equation  of  electrons  in  a  magnetic  field. 


(5b) 


Gas  Dynamic  Equations 

The  gas  dynamic  equations  consist  of  the  following: 

a.  The  differential  expression  of  the  mass  continuity  equation 

^u)  »  0  — ^  *  V  *  (AAu)  =  0  (for  variable  area  A) 


(6) 


2.  3fii.tzer,  L.,  Physics  of  Fully  Ionized  Gases 

3.  Allis,  W.  P.,  "Motion  of  Ions  and  Electrons,"  Handbuch  der  Physik,  Vol.  XXI, 
p.  397,  Eq.  15*3,  Springer-Verlag,  1956. 


b.  The  differential  equation  of  motion  which  is  essentially  a  statement 
of  Newton's  second  law  for  a  differential  volume  of  the  moving  medium 


-  ?p  +  JxB  +  +  /*g  +  P  y>V2u 


+  (1/3  )  (  V  •  u) 


(7) 


Energy  Equation 

The  differential  equation  for  the  generalized  energy  equation  for  a  com¬ 
pressible  fluid  flow  has  in  addition  to  the  internal  energy  term 


where  U  •  ^  C^dT^and  the  term  representing  the  work  done  by  motion 


1  t2 


heat  terms  per  unit  volume  due  to  Joule  heat  -  J  ,  due  to  heat  conduction 

2  (i  r  -i  2s 

k  v  T  and  due  to  viscositj  +  (  V  u)^ 


Thus 


^  [K +  -  *  *  + 2  *  *  + ? 


+  kV2T  +  \  f>  V£  Vu  +  (7u)# 


where  the  asterisk  denotes  the  transpose  of  the  tensor. 
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3.  MPD  Power  Generator  Steady-State  One-Dimensional  Flow  Equations 


By  applying  the  following  conditions  and  assumptions  to  the  generalized 
equations  governing  the  behavior  of  the  MPD  power  generator,  these  equations 
take  on  scalar  forms  which  can  be  readily  solved. 

a.  The  operation  of  the  generator  is  in  the  steady-state  (i.e.  =  0) 

b.  The  gas  flow  is  strictly  one  dimensional  which  means  that  all  variables, 
in  addition  to  the  gas  dynamical  variables,  vary  only  in  the  direction 
of  flow  (  x-direction) ,  but  not  over  the  cross  section  of  the  channel 

c.  The  applied  magnetic  field  is  in  the  z-direction  and  is  constant 

(i.e.  Bw  =  B) 
z 

d.  The  current  flow  is  possible  in  the  x  and  y  directions  only  (i.e.  =  0) 

e.  The  electrical  conductivity,  <j-q,  of  the  ionized  gas  is  constant 
(i.e.  o-Q  -  e( n^e+np/tp)) 

f.  The  electron  pressure  gradient  of  the  ionized  gas  is  zero,  i.e. 

’J.' 

g.  The  ionized  gas  is  electrically  neutral  (i.e.  n  =  n  ) 

e  •  p 

h.  The  effects  of  gravity,  viscosity  or  thermal  conductivity  of  the  gas  can 
be  neglected 

i.  The  ionized  gas  or  plasma  obeys  the  perfect  gas  law  (i.e.  p  =  ^RT) 

j.  The  ratio,  y,  of  specific  heat  at  constant  pressure  to  specific  heat 
at  constant  volume  of  the  gas  is  constant 

k.  The  ionized  gas  flow  through  the  generator  duct  is  free  from  shocks, 
i.e.  the  flow  is  free  from  any  regions  in  which  the  temperature^ pressure 
or  density  change  very  rapidly  with  distance 

l.  The  voltage  drops  at  the  electrode  sheaths  are  small  in  comparison  with 
the  generated  voltage  and  can  be  neglected. 
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Thus  the  generalized  equations  governing  the  behavior  of  the  MFD  power 
generator  can  be  reduced  to  the  following  scalar  forms: 

a.  Ohm's  Law  (neglecting  electron  pressure  gradients  and  ion  slip) 

1)  The  current  density,  Jx,  in  the  direction  of  flow  is  given  by 

J*"r75[e*  v(£y-u3)]  (5) 

2)  The  current  density,  J  ,  across  the  interelectrode  spacing  (perpen- 

y 

dicular  to  the  direction  of  flow  and  perpendicular  to  the  magnetic 
.  field) 

(10) 

where 

■/®*«t-yieB  (11) 

b.  Gasdynamic  Equations 
1.)  Continuity  Equation 


■ 0 

(12) 

puA  ■  O  u  A  ■  constant 
r  fo  o  0 

(13) 

2)  Equation  of  Motion 

*  JB 

r  dx  dx  y 

(Hi) 

Energy  Equation 

p  u  .  jj2  -  E  u  +  1_  j2 

r  dx  p  ax  o-o 

(15) 

where 

U  -  CyT  and 

J2  -  J2  +  J2 
x  y 

(16) 

(17) 
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d.  Equation  of  State 
Rn 

p  -  pRU  -  /5  -°T 


(18) 


where  Rq  ■  8.31I4.  Joules/°K  is  the  universal  gas  constant  and  Z  is  the 
atomic  weight  in  Kg(e.g.  Z  *  O.OOU  Kg  for  He). 

Differentiating  the  equation  of  state  (18)  with  respect  to  x  and  then 
substituting  into  the  energy  equation  (15)  we  get  the  following  expression 


P  u  !c  (CV  +  R)T  "  u  %  +  -ZT  j2  (19) 

0~o 

By  substituting  for  -jE  from  the  equation  o_  motion  (1U)  into  the  energy 
equation  (19)„and  since  =■  +  R  we  get,  after  rearranging  terms,  the  fol¬ 
lowing  expression  for  the  energy  equation 


uBJ  +  i-  J2 
y  O-o 


(uB  +  i-  J  )  J  + 
<T0  y  y 


-i-  J2 
0“o  x 


(20) 


A.  Solution  For  the  Constant  Velocity  Case 

In  this  case  u  ■  u  =  constant  and  the  MPD  power  generator  gas  dynamic 
0 

and  energy  equations  reduce  to  the  following  simple  form: 

a.  Continuity  Equation 

pk  *  p  QkQ  ■  constant  “  p  Q GqD0  *  y0GDQ  (21) 

b.  Equation  of  Motion 


(22) 
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c.  Energy  Equation 


I:  <V>  -  *  V0  VA 

1)  Segmented  Electrodes  Faraday  Generator 

In  this  case  ■  0  and  solution  for  the  electric  field  in  the  x 
direction  from  equation  (9)  yields  the  following  expression  for  ^ 


(23) 


£x  3  Pity-*** 


(2k) 


Substituting  this  value  of  £  into  the  expression  for  J  (Eq.  10) 

x  y 

we  get  the  following  expression 


Jy  =  ff*o(£y  “  uB) 


which  for  an  electric  load  factor 

en 

_  m  _  d  y  o 

*  ODen  circuit  voltage  "  uBD 

o 


e 


L  V  open 
oc 


can  be  written  as  follows 

■J  ■  -  0-quB(1  -  e^)  ■  constant  for  constant  0^ 


(25) 


(26) 


(27) 


Substituting  this  equation  (27)  for  J  in  the  equation  of  motion  (22) 

y 


we  get 


-u<r0B 


;(i- 


el) 


(28) 


Solution  of  equation  (28)  gives  for  the  pressure,  p,  at  any  distant  x 
the  following  expression 

P  *  P0  -  u<J-oB2(1  -  eL)x  (29) 

Now  from  the  equation  of  state  and  the  definition  of  Mach  number, 

2  2 

M  (viz.,  u  *  M  #RT),  we  can  derive  the  following  expression  for  the 
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density,  p  ,  in  terms  of  the  velocity,  u,  the  pressure,  p,  and  the  Mach 
number,  M,  namely  ... 

u 

By  substituting  into  the  energy  equation  (23)  the  expression  for  p 
given  by  equation  (30),  the  expression  for  J  given  by  equation  (27)  and 

o  y 


(30) 


the  relation  C  T 
P 


RT 


u 


-g —  and  then  rearranging  terms,  we 


y  - 1  (  y  -  dm‘ 

obtain  the  following  differential  equation  for  the  Mach  number,  M,  namely 


dM  .£U<roB  (l"eL)dx 


Po  -  u  <roB  (1  -  eL)x 


(31) 


Solution  of  this  differential  equation  gives  for  the  Mach  number,  M,  at 
any  distant,  x,  the  following  expression 


In 


s-  -  -f 


In 


Po  -  u<r0B  (1  -  eL)x 


In 


Po  -  u  cr  0B  (1  -  eL)x 


(32) 


where  M'  is  the  Mach  number  at  the  entrance  to  the  magnetic  field-current 
o 

interaction  region.  Thus  the  Mach  number  at  any  distance,  x,  is  given  by 


m  r  p0  i 

Mo  ’  LPo  “  U<roB"(l  - 


oy2 


P0  °/2 

(r) 


(33) 


where 


DL-  (  y  -  Dej/y 


(3b) 


The  gas  density,  p  ,  can  now  be  given  by  substituting  the  relation 
between  p  and  M  given  by  equation  (33)  into  equation  (30).  Hence, 


y  m  p  ot, 

— jr  ( — )  P 
u  p 


u  p 


a  -  i  p  i  -ot 

-  /V  r } 


m  2(1  -a)/a 

/Vtt) 


(35) 
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Now  from  the  continuity  equation  (21)  the  dimension  G  of  the  duct  must 
be  varied  according  to  the  following  equation 


M  2(1- <*)/<*  p  (l-oo  p  (a-i) 

°“Goij"Go(r)  mQ0  {r) 

O  **  *0 


e± 

p 


or 


0  .A,1"* 


m  2(1  -  a )/« 

(TT) 

o 


(36) 


The  power  output  is  given  by  the  following  integral 

jy£y  DGdX  (37) 

0 

Substituting  for  J^.  from  equation  (27),  for  £  y  the  expression  uBe^ , 
for  G  from  equation  (36)  and  for  dx  from  equation  (28)  we  get  the  following 

o  (flt"  1)  n 

ueLGoQoPo  C  p"}  jr) 

o  *o 

"  -  -^T  Po^o-'o  l1  -{vj\ 

"  -  Pq^Vo  \X  "  (  TT)  ]  (38) 


integral 


«e-/c 


Mow  the  generator  volume  is  given  by  the  integral 


L,  P,  -  («-l> 

Jwxx-  j  d0o0(  E-) 


* 

-u 


f 


(TftB£d  -  «T) 


v  _ 


] 


(3?) 


DoGoP< 


-JU-  f1  -  ( £-)a  — r 

1  L  p0  J  uo-  B2( 


(1  -  eT)e 


L'  L 
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Dividing  the  total  output  power  by  the  volume,  we  get  the  following 
expression  for  the  average  output  power  density 


we  “  -»2tf-0B?eL(1  "  eL} 


and  since  u2  -  M2  Jf  RT  and  T  -  — — I?-- — *- 
o  o  o  0  i  +  -f  -  l  m2 


,  hence. 


(Uo) 


*.  *  “.fA11  -  eL> 

,2 


-  s 


R 


M  ‘ 


So 


¥  R- 


TT32 I^v-1 

<:  O 

,2  <*-«“o 


(rA(l  -  eT) 


'  ('pT  7 5  Tso  ^oB  TTjEZ^t  eL(1  -  eL> 

2  o 


(Ul) 


or 


( y  -  dm, 


°  eT  (1  -  eT ) 


(Ms) 


The  inlet  power  is  given  by 

P  uA  C  Tc  -  uD  G  p  RT 
fo  o  p  So  o  of  o  ( 


Wo*1  *  “o  > 

Hence  the  power  conversion  efficiency  is  given  by 


(W) 


-W_  A  V  p  ' 


M  2 


T«rTT3^--m^- 

since  from  equation  (33)  we  have 

p  u<TB2(1  -  eT)L  M  2/* 

p  p  v  ?r; 

*o  *o  r 


(U3) 


(UU) 


MPD  Solution 
Constant  Velocity  Case 


Ratio  of  Mach  Number  versus  Load  Factor 


MPD  Solution 
Constant  Velocity  Case 


Power  Density  vs  Inlet  Kach  Number 
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Fffirlency  vs  Inlet  Mach  Number 


-  3P 


where  and  are  the  pressure  and  Mach,  number  at  the  end  of  the  magnetic 
field  -  current  interaction  region  and  L  is  the  length  of  the  dftfct  cor¬ 
responding  to  exit  Mach  number,  Mr,  and  is  given,  by 


L  a 


r  m  2/*~ 

p 

D 

1  -  (~) 

V 

.  0 

u  r0Bz(i  -  ,p 

1  -2 
.  po. 

(W>. 


In  the  following  figures  the  solution,  equations  presented  above  for, the 

constant  velocity  segmented  electrodes  Faraday  generator  case  are  presented 

in  universal  dimensionless  form.  Thus,  for  example:,  if  we  fix  the  electric 

load  condition  (viz.,  eT)  and  the  exit  to  inlet  area  ratio  (viz.,  Gy/G^)  we 

determine  the  ratio  of  inlet  to  exit  Mach  number  (vix. ,  M  /M  ),  the  ratio  of 

o  r 

exit  to  inlet  pressure  (viz.,  Pr/pQ),  the  ratio  of  exit  to  inlet  temperature 
(viz.,  Ti/Tq)  and  the  generator  length  parameter  (viz.,  Lu  (TB/p^)..  Next 
by  choosing,  the  inlet  Mach  number,  Mq,  temperature,  T  ,  and  pressure  pQ  we 
can  now  calculate  the  generator  length,  L,  generator  efficiency,  ,  and  the 
generator  average  power  density,  wg,  (in  terms  of  the  power  density  factor 

Vs,";®2)  ho •  V1  *  mo2)- 

B.  Constant  Area  Duct  Segmented  Electrodes  Faraday  Generator  with  Constant 
Electrical.  Conductivity,  Q~o,  and  Constant  Internal  Voltage  Drop  Between 
Each*  Pair  of  Electrodes  Along  the  Length  of  the  Dust  ” 


In  this  case  the  current  density  in  the  x  direction,  J  is  zero  and,  hence, 
the  electric  field  in  the  flow  direction  is  given  by  equation  (2U),  viz., 

ex-  js  ( -  «b>  (a> 

And  the  current  density  in  the  y  direction  is  again  given  by  equation  (25),  viz.. 


Jy  ■ 


(25) 
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If  the  condition  is  further  observed  that  in  the  case  of  non-equilibrium 
ionization  conductivity  resulting  from  electron  heating*  it  is  necessary  to 
maintain  the  internal  electric  field,  and  hence  also  the  current  density, 

Jy,  constants  along  the  duct  length, then  Ohm's  law  (or  equation  2$)  becomes 

=  -  Cg-  (uB  -  gy)  =  -<ro  £e  =  -  <T0%B(1  -  eL)  =  constant  (16) 

For  the  constant  area  duct  and  under  the  condition  of  Ohm’s  law  in  equation 
(16)  the  continuity  equation,  the  equation  of  motion,  the  energy  equations  and 
the  equation  of  state  now  take  on  the  following  forms,  respectively* 


a.  Continuity  Equation 

p  u  B  p ouq  =  constant 

b.  Equation  of  Motion 

/0u  5  +  4  *  J  B  =  -  u  (h  B2(l  -  eT)  =  constant 
I  o  o  ax  dx  y  o  o  K  V 

c  .  Energy  Equation 


(17) 


(18) 


/Vo  cEc 


■cy^Dsr +  ? 


u2  =  (uB  +  i-  J  )  J  =  £  J 


(To  y  7  7  7 


<Vc 


uB  -  uqB(1  -  eL)  j  -  eTv 


u  -  uo(l  -  eL) 


-i  j 

)  uo(ToB2(1  -  eL)  (19) 


d.  Equation  of  State 


.2 


p  *  z>RT  ■  P ■  --A1  =  ft  u  (since  u2  •  M2yRT) 


m 


*  See  paper  by  M.  E.  Talaat  "Magnetoplasmadyriamie  Electrical  Power  Generation 
with, Non-Equilibrium  Ionization,"  presented  September  6,  1962  at  the  Symposium 
•  on  Maghetoplasmadynamic  Electrical  Power  Generation  at  Kings  College,  Newcastle 
Upon  Tyne,  England.  Paper  under  publication  in  Advanced  Energy  Conversion  Journal, 
1963j  see .  also  Martin  Semi-annual  Report  MND-2939. 
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By  substituting  the  expression  for  the  pressure,  p,  from  equation  (50) 
into  equation  (48)  we  get 


fin  $- 

r  o  o  dx 


-  -  ".n^1  -  ei> 

Now  by  multiplying  equation  (51)  by  u  -  uq(1  -  e^)^J 


(1  -  eT)  and  equating  the 


lefthand  side  of  the  result  with  the  left  hand:  side  of  equation  (U9)  we  get 
equation  (52)  which  gives  the  differential  relation  between  the  velocity,  u, 
and  the  Mach  number,  M,  viz., 

t  -  "»(1  -  *  [(1  *  y?»”]  ■  ®  Ifety  * .  *}  ”2]  m 

V  RT  T  t  * 

Now,  let  uq(1  -  e^)  =  ue  and  substitute  g  -g-  ?or  ***  equation  ("5  2')>  hence 

r 

(u  -  ue)d(u  +  f£)  =  d  l-y^-RT  +  £  u*  =  d(CT  +  \  u2)  =  CpdT  +  udu  (53) 


(u  -  u  )  du  +  SdT  -  ^  du  I  »  C  dT  +  udu 

e  L  u  /.  u  j  p 

and  since  Gp  =  +  prearranging  terms;  iii  equation  (54)  we  get 


u  )  + 
e 


2 

u  u  du  +  u 
e 


u  Cb  +  uu  R  dT  =  0 
v  e 


Now  let'  T  =  vu  hence, 

dT  =  vdu  +  udv  ■ 


Substituting  from  equations  (56)  and  (57)  into  equation  (55)  and  eliminating 
T  and  dT  from  equation  (55)  we  get  after  rearranging  the  terms 

(Cpv  +  ue)du  ■  -  (uCy.  +  ugR)dv  (58) 


1  °pdv _ 

C  C  v  +  u 
P  P  e 


!  cVdu 

^(^u+'ftue 
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The  relation  between  the  velocity,  u,  and  the  temperature,  T,  or  Mach 

number,  M,  can  be  readily  obtained  by  integrating  equation  (5>9)  subject  to  the 

boundary  conditions  that  at  x  =  0,  v  -  v  ,  u  =  u  ,  T  =  T  and  M  =  M  .  Thus, 

o  o’  o  o  * 


C  v  +  u 

J2 _ £_ 

C  v  +  u 
pi  o  e 


Fc.-u 

+  Ru  1 

1  V  0 

e 

C..U  +  Ru 

l 

6  J 

t  y 

and  since  v  »  -  ,  Cl  =  -  rt  = 

u  *  p  y  -  l 


(  i  -  1)M 


2  +  ue  = 


( *  -  i)  m; 


2  +  ue 


T  u 

,  hence,  C  v  =  C  -  =  rtr—rvo 

*  *  p  p  u  cy-ijM^ 


CTru  +  Ru 
Vo  e 


CT.u  +  Ru 
V  .  e 


1 1  T5  . 

let  the  velocity  ratio  s  =  —  ,  and  since  ^-  ■  (  6  -  1)  and  recalling  that 

o  V 

u  =  (1  -  e,)u  we  get  by  substitution  into  equation  (61)  the  following  relation 

6  Jj  O 

between  the  velocity  ratio  and  the  Mach  numbers 


Note  that  at  eT 


rs  i  i  m  + ( # -uu  -  eLn* 

Lm2  LJ  J_M2  JLs  +  (  X  -1)(1  “  eL)J 

s  i  x  2/c  y  +  i) 

Note  that  at  eT  =  l."—*  =  — *  ~~s  or  s  =  (  -*-)  which  checks  the 

L  M2  M2  s*  Mo 

o 

isentropic  flow  relation  as  it  should. 

Next,  solution  of  the  equation  of  motion  (equation  51)  gives  the  relation 

between  the  axial  distance,  x,  measured  from  the  inlet  to  the  magnetic  field- 

current  interaction  region  and  the  corresponding  velocity  ratio,  s  »  and  Mach 

number,  M,  for  a  given  inlet  Mach  number  and  a  given  electric  load  factor,  e^, 

thus  by  integrating  equation  (63) 

_  -  -  _  .2 


d  f.  ,  s  1  (To® 


(1  -  eL) 
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subject  to  the  boundary  conditions  that  at  x  =  0  s  ■  1  and  M  »  M  we  get 


C1  +  -  s  (1  +  (i  _  eT )  x 


If  the  Mach  number  at  the  exit  is  M  and  the  corresponding  velocity  ratio 
is  er  the  generator  active  length,  L,  can  be  obtained  from  equation  (61)),  thus 


LOqB  „ 

f>0\ 


i  r(i +jh 

“  eL^  [_  ^  M' 


O  v  V 


In  equation  (65)  once  Mp  is.  chosen  the  corresponding  s^  can  be  obtained  from 
equation  (62)  which  gives  the  relation  between  M  and  s  for  a  given  inlet  Mach 
number,  Mq,  and  a  given  electric  load  factor,  eL.  Equation  (65)  can  thus  be 
used  to  determine  L  for  a  given  set  of  operating  conditions  (Mq,  e~,  Tq,  Mr  and 
8r)  or  determine.  the  exit  Mach  number  and  velocity  ratio  once  L  is  chosen. 

The  output  power  is  given  by  the  following  expression 


L 

■  /  “VSy* 


Substituting  for  J^. Sydx  from  equation  (U9)  into  equation  (66)  the  expression 


for  the  output  power  reduces  to  the  following 


W  =  DG 
e 


r  2 

v°  l , 4 


u  ,M 

o*  o 


r  2  2  1 

V>0U°U°  f  df(^l)M2  +  2J 

> lo 

SA"o  "o  [tTTTJ  <7-3>  *1  -  »] 
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The  expression  for  average  output  power  density,  wg>  is  next  obtained- by 
dividing  the  expression  for  output  power } equation  (67),  by  the:  volume  DGL  and 
using  the  expression  for  L  as  given  by  equation  (6£) .  Thus, 


the  expression  for  the  average  output  power  density  can  be  written  as  follows 


(70) 


Thus  for  a.  given  gas  at  a  given  stagnation  temperature,  T  ,  inlet  and  exit 
Mach  numbers,  Mq  and  M  ,  electric  load  factor,,  e. ,  and  exit  to  inlet  velocity, 

sp,  (see  equation  (62)  for  dependence  of  sr  on  M  ,  Mq  and  e^)  we  can  evaluate  the 
average  power  density,  we> 

Finally  the  expression  for  power  extraction  efficiency  is  given  by  the  ratio 
of  output  power  to  total  enthalpy  at  the  duct  inlet  (see  case  A  and  Ref.  1) 


n 


i  + 1 

(  X  -  i)  ? 


( *  -  Dm? 


(71) 
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where  again  the  exit  to  inlet  velocity  ratio,  sr,  is  determined  from  equation 
(62)  once  My,  Mq  and  e^  are  specified. 

In  the  following  figures  the  solution  equations  presented  above  for  the 
constant  area  duct  segmented  electrode  Faraday  generator  case  are  presented 
in  universal  dimensionless  forms.  The  exit  to  inlet  velocity  ratio  is  thus 
plotted  versus  the  inlet  Mach  number,  Mq,  for  a  constant  exit  Mach  number 
My  *  1  (the  choke  condition  since  in  this  case  the  flow  is  choke  limited)  and 
electric  load  factor  e^  of  0.5,  0.7  and  0.9.  Similar  figures  are  given  for  ; 

p  p  p 

the  exit  to  inlet  temperature  ratio',  (Tr/T())  =  s£(M^/Mp  and  the  exit  to  inlet 

?  p 

pressure  ratio,  (Pj/Pq)  -  s,r(M^Ar)  again  for  M^  *  1  and  eL  -  0.5,  0.7  and  0.9. 

2 

The  length  parameter  (L  <TqB  / n  uq)  is  plotted  versus  the  inlet  Mach 

number,  Mq,  for  an. exit  Mach  number,  -  1,  and  load  factors  eL  ■  0.5,  0.7 

and  0.9.  Similarly,  the  energy  conversion  efficiency,  Tj  ,  is  plotted  versus 

the  inlet  Mach  number,  M  ,  for  an  exit  Mach  number,  M  ■  1,  and  load  factors 

o  *  p  ' 

p 

*  0.5,  0.7  and  0.9.  The  power  density  parameter  ( We/Cp Tsq <TqB‘: )  is  also 

plotted  ■  versus  the  inlet  Mach  number,  M  ,  for  an  exit  Mach  number,  M  *  1  arid 

0  r 

load  factors  e^  ■  0.5,'  0.7,  0.8,  0.85  and  0.9.  Note  that  maximum  (100$) 
efficiency  for  any  load  factor,  e^,  is  approached  as  the  inlet  Mach  number, 

Mq,  approaches  zerofwhile  zero  efficiency  is  approached  as  M  approaches  Mach 

l 

one.When  the  inlet  Mach  number,  Mq,  approaches  CO  the  efficiency, ,  for  any 
load  factor,  e^,  approaches  the  value  given  by 

1  -  h  sL  (72) 

where  is  the  exit  to  inlet  velocity  ratio  for  M O'1*‘co  and  is  given  by 

,  X  0.375  r  )  0.625 

s*  -  {  (  t  -  D(1  -  eL)|  (1  +  (  -  1)(1  -  eL)j  -  (*-1X1  -  eL) 

(73) 
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Thus  for  a  monatomic  gas  at  load  factors  e^  =>  0.5*  0.7,  0.8,  0.8$  and  0.9  the 
exit  to  inlet  velocity  ratios  s  ^  are  0.14,585,  0.U13,  0.3733,  0.3U8  and  0.3113, 
respectively  and  the  efficiencies,  ,  approached  as  Mq  approaches  00  are 
0.160,  0.318,  0.UUU,  0.5185  and  0.613,  respectively. 

p 

The  power  density  parameter  (w^CpT^tJ^B  ),  on  the  other  hand,  approaches 
zero  as 'the  inlet  Mach  number  approaches  zero.  As. the  inlet  Mach  number  approaches 
one  the  power  density  parameter  (we/GpTgo  (T^B  ). approaches  the  following  expression 


C  T 
p  so 


M  -*>1 


(1  -  eL) 


+  (*-l)(l-eL) 


} 


(7U) 


o 

Thus  for  a  monatomic  gas  at  load  factors  e^  =  0.5,  0.7,  0.8,  0.85  and  0.9  the 
power  density  factors  approached  as  Mq  approaches  unity  are  0.125,  0.105,  0.08, 
0.06375  and  0.0U5,  respectively.  As  the  inlet  Mach  number  approaches  infinity 
the  power  density  factor  approaches  the  limit 


M-^oo 


C  T  (T’B' 
p  so  o 


U  -  eL) 


y  (y-D  -  ( s  +  Dsi 
_)  _____ 


o 

Thus  for  a  monatomic  gas  at  load  factors  e  ■  0.5,  0.7,  0.8,  0.85  and  0.9  the 

L 

power  density  factors  approached  as  Mq  approaches  infinity  are  0.300,  0.281, 
0.221,  0.175  and  0.122,  respectively. 


(75) 
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G.  Constant  Area  Duct  Segmented  Electrodes  Faraday  Generator  with  Constant 
Electrical  Conductivity,  (T'm  and  Constant  Electrical  Load  Resistance  for 
Each  Electrode  Pair  Along  the  Duct 


In  this  case  the  current  density  in  the  x-direction  is  zero  and  the  current 
density  in  the  y-direction  is  again  given  by 

Jy  -  <r^(  £y  -  UB)  -  -  0-0uB(l  -  ^5L)  (76) 

Now  if  we  let  (  S^/uB)  -  e^  *  constant  (77) 

for  each  pair  of  opposite  electrodes,  then  the  expression  for  the  current  J  can 
be  written  as  follows 

■  fjftt  'eL>  (78> 

which  when  solved  for  (  6y/Jy)  yields  the  following  expression 

(6/y  -  -  eL/(To(l  -  eL)  (79) 


Equation  (79)  in  essence  states  that  the  load  condition  expressed  by  equation 
(77)  exists  when  the  load  resistance  is  maintained  constant  for  each  pair  of 
electrodes  along  the  generator  duct. 

The  equations  governing  the  behavior  of  the  MPD  generator,  under  these 
conditions,  then  taka  the  following  forms 

1)  Ohm's  Law 

Jy  -  -  0^uB(l  -  eL)  (80) 

2)  Continuity 

p u  ■  “  constant  (81) 

3)  Equation  of  Motion 

f> 0%  « *  %  m  V5  "  “  * ^o®2*1  "  eL)  (82) 
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ii)  Energy  Equation 


1  i  - H - +  i-  »  u' 

“ o  o  ax  ( %  -i)>r  * 


2  -  eL) 


(«3) 


5)  Equation  of  State 


(89) 


Solution  of  these  equations  yields  the  following  expression  for  the  velocity 
ratio,  s,  in  terms  of  the  load  factor,  e^,  the  inlet  Mach  number,  Mq,  and  the 
Mach  number,  M,  at  any  point  at  distance  x  from  the  inlet  to  the  magnetic  field  - 
current  interaction  region 

»2/  .,,2  .  i_  \  (c/b) 


rwK>b)  i 

Uo  (aM^  +  b)  _ 


where 


a-  X(y-D(l-eL) 
b  -  2  -  (If  -  l)eL 

c  •  t  -  (%-  l)eT 


(86) 

(87) 

(88) 


Solution  for  the  relation  between  the  distance  x,  the  load  factor,  e^,  the 
inlet  Mach  number,  MQ,  the  Mach  number,  M,  at  the  distance  yields 

_  1  f  1  1>  Odff!)  /nMWj_b) 

“  elP  m!  1?  b2(l-eT)  M2(aM^  ♦  b) 


f>o\ 


(89) 


The  expressions  for  the  output  power  and  efficiency^  are  given  in  the 
following  and  are  derived  by  following  similar  arguments  as  those  used  in 


deriving  parallel  expressions  for  cases  A  and  B 


Thus 


-  1*0  - 


»e  ■  -®  /■>,£,  <“■ -  ®/Wo  J  d  (  — i_,  ♦ 
o  i,mq  L  (  (y-i)M2 

■  'tv  "o  [[77V  *  *  i  -  °] 


V  * 


1  +I 


( W>°  7H)  [7-  *  ^1 1  -^1 7\  '  •*] 

<®/’ouo>  x^rOTo  [x  *  ^  Mo]  7 


®/>o“o  SM 


where  the  efficiency  ^  is  again  given  by 

? 


^  T 


(90) 


(91) 


The  average  output  power  density,  wg,  is  obtained  by  using  equation  (90)  and 
(89)  and  is  given  by  the  following  exprescion 


?  b(1  ~  eL} _  (92) 

-  -j)  -  (  if  +  1)  A  n  b 

where  V)  is  given  by  equation  (91)  and  s  by  equation  (85). 


VsonB 
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D.  Constant  Area  Duct  Segmented  Electrodes  Faraday  Generator  with  Constant 
Electrical  Conductivity,  and  Constant  Output  Voltage  Along,' the  Pact 


In  this  case  the  load  voltage  is  kept  constant  between  each  pair  of  elec¬ 
trodes  along  the  duct.  Thus,  £  ■  constant  =  eLu0®0  an<^>  hence, 

Jy  "  -W(S  "  eL}  (93) 

Using  this  expression  for  the  current  density  in  the  y-direction  the  equation 
of  motion  and  the  energy  equation  take  on  the  following  forms 


1)  Equation  of  Motion 

p  ouo  It (u  +  j^r)  m  Vz "  -  (Touob2(s  *  eL) 


(9h) 


2)  Energy  Equation 


PoUok  {/yW  +  2  I  =  £yJ3 


-uo  -  eL> 


(95) 


Equations  (9b)  and  (S>5)  are  the  same  as  equations  ( 20)  and  (21),  respectively, 
in  Ref.  1  and  the  solution  of  this  case  has  been  presented  in  that  reference. 
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